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ABSTRACT 

Black holes release energy via the production of photons in their accretion discs but 
also via the acceleration of jets. Wc investigate the relative importance of these two 
paths over cosmic time by determining the mechanical luminosity function (LF) of 
radio sources and by comparing it to a previous determination of the bolometric LF 
of active galactic nuclei (AGN) from X-ray, optical and infrared observations. The 
mechanical LF of radio sources is computed in two steps: the determination of the 
mechanical luminosity as a function of the radio luminosity and its convolution with 
the radio LF of radio sources. Even with the large uncertainty deriving from the 
former, we can conclude that the contribution of jets is unlikely to be much larger 
than 10% of the AGN energy budget at any cosmic epoch. 

Key words: black hole physics — galaxies: active — galaxies: jets — 



1 INTRODUCTION 

Matter can accrete onto a black hole (BH) only if it releases 
a fraction e ~ 0.06 — 0.4 of its rest-mas s energy, where e 
depends on the BH spin ([Bardeenl Il970l ). In the standard 
IShakura fc SvunvaevI (119731 ) model, the energy is dissipated 
by viscous torques in the accretion disc and radiated. Ac- 
cretion from a luminous disc pro vides a physical mo del to 
explain the luminosity of quasars (|Lvnden-Belllll969l ). 

Several authors have computed the total energy 
radiated by BHs over cosmic time an d have com- 
pared it to the local B H mass density dSoltanl 1 19821: 
Chokshi fc Turner 1 992: Yu fc Tremainell20Q2l : iBarger etal] 



20051 : Hopkins et al.l 120071 '). The two are in good agreement 



if matter is turned into light with a canonical efficiency of 
e ~ 0.1. This has been used as an argument to infer that 
most of the BH mass in the Universe was accreted lumi- 
nously, but it only proves that > 25% of the BH mass in the 
Universe was accreted luminously, since e could be as large 
as e ~ 0.4 if most BHs were maximally rotating. 

In fact, accreting BHs ('active galactic nuclei', AGN) 
produce not only light but also jets of matter, which are 
radio-luminous because of the synchrotron radiation from 
ultrarelativistic electrons accelerated in shocks. Excluding 
objects beamed toward the line of sight, only a small fraction 
of the luminosity that is radiated by an AGN comes from 
synchrotron radiation. However, the synchrotron power rep- 
resents only a small fraction of the jet mechanical luminosity, 
most of which may be used to do work on the surrounding 
gas. 

Moreover, at low accretion rates (M. ^ O.OIA/Edd, 



where Z/Edd = eMEddC^ is the Eddington luminosity), the 
accretion disc is not dense enough to radiate efficiently; the 
disc puffs up, and the energy that needs to be removed 
to allow the accretion may be carried out more easily by 
jets. Although the physics of this picture are still specu- 
lative, AGN that channel a large fraction of the accretion 
power into jets while sho wing little emission from an accre- 
tion d isc are observed (eg lDi Matteo et al.ll2003l : lAllen et al.l 
l2006h . These radio sources are less powerful than quasars 
but more common due to their longer duty cycle. For ex- 
ample, around two-thirds of brightest cluste r galaxies are 
radio galaxies (jBurnsI [l990l : iBest et al.ll2007l ). In contrast, 
only one g alaxy in 10^ contain s a quasar {Mb < —23) 
at 2 ~ (|Wisotzki et al.l 1200 ih . Finally, mechanical en- 
ergy is thermalized in the intracluster medium more ef- 
ficiently than luminous energy is. The observational evi- 
dence that the mechanical heating by AGN is important to 
solve the cooling-fl ow problern in gal a xy groups and clusters 
is getting strong jBest et al.l [ioOSal: iDunn fc Fabianll2006l : 
iRaffertv et~al]|200g : lMagliocchetti fc Briiggenll2007l : see also 
Cattaneo et al. 2009, and references therein). 

For these reasons, it is important to compare the me- 
chanical and radiative output of AGN. This is the goal of 
the curre nt paper. This issue has also recently been ad- 
dressed bvlShankar et aD (|2008ah . lK5rding et al.l (j2008l l and 
iMerloni fc Hein3 1 20081 ) . We adopt a different approach to 



these authors and produce results that are qualitatively sim- 
ilar. The layout of our paper is as follows. In Section 2, we 
analyse how we can use radio data to infer the jet mechanical 
power. We shall see that two different approaches give differ- 
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ent Pjet(iradio) relations. We consider both, and use the dif- 
ference between the results obtained from the two relations 
to provide an estimate of the uncertainty. We convolve these 
relations with the radio LF of radio sources, (?!)(Lradio), to es- 
timate the mechanical LF, 0(Lradio(fjet)), first in the local 
Universe (Section 3), then at different redshifts (Section 4). 
In both cases, we integrate over luminosity to determine the 
mechanical power per unit volume, which we compare with 
the radiative power per unit volume from luminous AGN. 
In Section 5 we summarize the results of comparing the me- 
chanical LF of AGN t o the bolometr i c lum inosity function 
of AGN determined bv lHopkins et all (|2007l ) and we discuss 
the implications of our results. 



2 THE MECHANICAL LUMINOSITY OF A 
RADIO SOURCE 

Obtaining an estimate of the mechanical power of a ra- 
dio source is an inherently difficult problem. The observed 
monochromatic radio luminosity measures only the fraction 
of the jet power that is currently being converted into radia- 
tion. That frac tion is small (typically between 0.1 and 1%; cf. 
lBicknelllll995l ) and changes during the lifetime of the radio 
source, since the radio luminosity of a growing radio source 
first increases and then drops as the source expands into 
a pro gressively lower density environment (eg lKaiser et al.l 
1 19971 ). Nevertheless, it is reasonable to expect that radio and 
mechanical luminosities should show at least a broad degree 
of correlation on a population basis. 

Estimates of the mechanical po wer of radio sources have 
followed two approaches. The first (jWillott et al.|[l999l ) uses 
the minimum energy density Wmin that the plasma in the 
radio lobes must hav e in order t o emit the observed syn- 
chrotron radiation reg |MiIevl[l980l ). With this approach, the 
jet mechanical power is Z/moch ~ UminV/t, where V is the 
volume filled by the radio lobes and the radio-source life- 
time t is given by the ratio between the jets' length and 
the hotspots' advancement speed. The largest sources of 
uncertainties are; (i) the nature of the jet plasma (electron- 
positron or electron-proton?): the value of Mmin is larger if 
the lobes contain a hadronic component in addition to the 
synchrotron radiating particles (relativistic electrons and/or 
positrons); (ii) the lack of observational constraints on the 
low-frequency cut-off of the electron energy distribution: for 
a synchrotrom spectrum oc with radio spectral index 
a < —0.5, limin is larger when the lower cut-off frequency 
takes a lower value and thus there is mo re energy in the syn- 
chrotron spectrum. IWillott et all (|l999l ) derive the relation 



1.4 X 10^ 



Li. 



4 GHz 



3 X 10"»/^ 



/2( 

.T I 



Ll5 



Vl028WHz~isr-i 



6/7 



w, 



(1) 



wher e fw '•-^1 — 20 inc o rpora tes all the unknown fac- 
tors. iBlundell fc Rawhngsl (|200GI 1 argue for fw — 10 for 



Fanaroff fc Rilev (Il974f) class II sources (FRIIs), while 
Hardcastle et al I (|2007l ') suggest /w = 10 - 20 for FRIs. 
We convert the luminosity at 151MHz, Lisimhz, into a lu- 
minosity at 1.4 GHz, Li.4GHz, where the local radio LF is 
best-determined. This will also allow us later to compare 
Eq. ll]) with another determination of Lmcch by Best et al. 
(2006, 2007). For the conversion we assume a spectral index 
Q = -0.8. Using fw = 10, Eq. (Q gives: 



, 1025 W Hz- 
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Notice that Lisimhz is given in WHz ^sr ^ while Li.4ghz 
is given in W Hz ~ ^ to r espe ct the differ e nt con ventions used 
bv lWillott et al.l (|l999l ') and lBest et"all (120061 ). so there is a 
factor of 47r entering the conversion. 

A second approach is to infer Lmcch from the mechani- 
cal work that the lobes do on the surrounding hot gas. The 
expanding lobes of relativistic synchrotron-emitting plasma 
open cavities in the ambient thermal X-ray emitting plasma, 
which advances in X-ray imaging capabilities now allow to 
be imaged in detail. The minimum work in inflating these 
cavities is done for reversible (quasi-static) inflation and 
equals pV , where p is the pressure of the ambient gas. 
iBest et ahl (|2006l ) derived a relation between radio and me- 
chanical luminosity based upon this estimate for the en- 
ergy associated with these cavities, combined with an es- 
timate of the cavit y ages from the buoyancy timescale (from 
iBirzan et al.ll2003 ). Comparing the mechanical luminosities 
of 19 nearby radio sources that have associated X-ray cav- 
ities with their 1.4 GHz monochromatic radio luminosities 
leads to a relation 

imoch = (3.0 ± 0.2) X 10=''' / f — ^ 



Vl025WHz-iy ^'^^^ 
broad ly in agreement with that derived by iBirzan et al 



(2004). In Eq. (|3]), the factor /, incorporated by Best et al 



((20071 ) ■ accounts for any systematic error in estimating the 
mechanical luminosities of the cavities. In particular, pV is 
likely to be an underestimate of the energy needed to inflate 
a cavity: the enthalpy of the cavity is ^^pV — 4pV for the 
relativistic plasma in the radio lobes, suggesting that / ~ 4 
may be appropriate. Some authors ha ve even argued for 
mechanical energie s in excess of lOpV^ (|Nusser et al. I [20061 : 
iBinnev et al.|[2007[ ) due to additional heating directly from 
the jets. For / = 4, Eq. (|3j gives: 



= 1.2 X 10^ 
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Eqs. [2] and |4] are in excellent agreement at Li.4ghz ~ 
lO'^^WHz^^, but the relation derived from the minimum- 
energy argument has a steeper radio luminosity dependence 
than the relation derived from the X-ray cavities. Neverthe- 
less, given the totally independent approaches used to derive 
Eqs. [2]and|4l and the uncertainty factors in both equations, 
the degree of consistency is encouraging. Indeed, the dis- 
crepency between the two relations may simply reflect the 
fact that Eq. [5] is determined from powerful, radiatively ef- 
flcient radio sources (mostly I/i.4ghz > lO^^WHz-i) while 
Eq. |4]is derived predominantly from lower luminosity, radia- 
tively inefficient sources. The mean relation between Lmoch 
and Lrad may well be different in these two different regimes, 
in which case a combination of Eq.[5]at high luminosities and 
Eq. |3]at low luminosities would be most appropriate. 

A third approach towards estimating mechanical 
luminosities of radio sources has been developed by 
iMerloni fc Heinj (|2008t ). These authors use the de-beamed 
radio core emission as a measure of the jet kinetic lumi- 
nosity, based upon the analogy with X-ray binary sources 
and the so-called Fun damental Plane relation for black holes 
(|Merloni et al.ll2003l ). This approach uses the radio LF of 
flat-spectrum (ie core-dominated) radio sources as a measure 
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of the LF of radio cores. It then requires assumptions about 
the statistical de-beaming of radio sources (ie the distribu- 
tion of Lorentz factors of jets), and about how to correct for 
the radio cores that are missed from the flat-spectrum radio 
LF because their radio sources are dominated by extended 
steep spectrum emission. These factors can be reasonably 
estimated for moderate to high radio luminosity sources in 
the local Universe, but they are not so well constrained at 
low radio luminosities, where the bulk of the jet mechanical 
power is produced, nor at higher redshifts. We therefore do 
not consider this approach here, b ut we do compare our re - 
sults with the results obtained by iMerloni fc Heinj (|2008l ) 
in Section [5l 



3 THE LOCAL MECHANICAL LF 

In order to derive a mechanical LF for radio sources in 
the nearby Universe, we must convolve Eqs. [2] and |4] with 
the lo cal radio LF. Here we adopt the local 1.4 GHz radio 
LF of lBest et all (|2005al ). which is derived from the Sloan 
Digital Sky Survey spectroscopic sample and is fully con- 
sistent with other recent determinations of the local radi o 
LF (eg iMachalski fc Godlowskil I2OO0I : ISadler et al. I l2002t ). 
Throughout this paper we follow the convention of defin- 
ing the LF as the number of objects per unit vol- 
ume and logarithmic-luminosity interval (eg iHopkins et al.l 
I2OO7I ). so the number of sources with luminosity between 
L and L + dL is ((f)/lnlO)dL/L. With this definition, the 
iBest et al.1 (12005^ 1 LF can be be parameterised using the 
double power-law model 



0(L) = 0. 



l:) +(1: 



(5) 



where L = I/i.4ghz (we wrote L without any subscripts 
because we shall use Eq. [5]to model other LFs). The best- 
fit parameters are 0, = 10"''-^Mpc-^, L. = lO^'^'^'^ W Hz-\ 
OL = 0.57 and /3 = 2.31. Combining Eq. ^ with Eqs. [2] 
and[4]gives the blue and the red dotted curves in the 2 = 0.1 
diagram of Fig. 1, respectively. 

The mechanical power that is released per unit volume 

is 



Pmcch 



in 



4){L) dL 
InlO T 



of radio-loud AGN would exceed that of galaxies (or super- 
massive black holes) capable of hosting them. For example, 
if the slope of the radio LF were to remain unaltered down 
to Li,4 GHz = lO^'^WHz-i then the local space density of ra- 
dio galaxies would exceed that of galaxies integrated down 
to Mb ~ -15. 

Considering that only massive galaxies with massive 
black holes have a significant probability of harbouring a 
radio source powered by an AGN, the local space den- 
sity of radio sources matches that of supermassive black 
holes with M, > 10^ M r (ie ~ 5 x 10"^ Mpc~^ according to 
IShankar et"ai]|2008bl ) if the radio LF is extrapolated down 
to 1/1.4 GHz ~ lO'^^'^W Hz~^. This sets a strong lower limit 
on the luminosity to which the radio LF can be extrapolated, 
and so for a conservative calculation we evaluate the inte- 
gral in Eq. © by adopting lO^^'^W Hz" V^* - 10"^ as the 
lower extreme of the integration interval. With this choice, 
the mechanical power per unit volume from Eq. Q, PmochJ 
is about ten times larger than the mechanical power per unit 



volume from Eq. ((2)1 , p! 



synchro 
mcch 



The discrepency would reduce 



to a factor of three if a limit of lO^^WHz"^ were adopted 
instead. 

We want to compare these values to the radiative power 
per unit volume from luminous AGN. iHopkins etal] (|2007l ) 
made the first attempt at determining the bolometric LF 
of AGN by combining hard X-ray, soft X-ray, optical and 
infrared data. They found that the double power-law model 
in Eq. ([SJ fits their results with 0,, L*, a and /3 dependent 
on redshift (here L is the bolometric luminosity -Lboi, not 
the radio luminosity). The thick solid lines in Fig. 1 show 
their best fit at different z. We use their best fit to the 
bolometric LF at z = 0.1 to estimate the power per unit 
volume radiated by AGN in the local Universe: 



Prad = 



(jjjLhoi) dLbol 
InlO Lboi 



(7) 



(6) 



We find that pj^"c'if'^°//^rad — 0.018 (the level of the dashed 
line in Fig. 2). The value of Pmcch/Prad depends on the ra- 
dio luminosity Lmin that is taken for the lower bound of 
the integration interval in Eq. ((6} (solid line in Fig. 2). For 
Lmin/i, ^ W~^, PmJch/Prad ^ 0.2. For p^^^h to be compa- 
rable to Prad, the power law of the radio LF would have to 
extrapolate down to //i.4ghz ~ lO^^WHz"^ (Fig. 2), which 
is well below the allowed limits. 



= A 



Vl025 W/ 



dx 



where ^ = 4.8 x IO^Lq and rj = 0.85 for Eq. Q, and 
A = 1.0 X 10"Lq and rj = 0.4 for Eq. dH). If L^cch is 
computed with Eq. then p^^^t"" = 1-6 x lO^L© Mpc"^ 
{ILq = 3.9x 10^*^ W). If Lmoch is computed with Eq. Q, the 
problem is more complicated because in that case the inte- 
gral in Eq. ((6)1 diverges at 2: = 0. This is because the double 
power-law model in Eq. ((5} cannot be extrapolated down to 
L —r 0, as can be easily seen: the faint-end slope of the local 
radio LF is steeper than that of bo th the low luminosity end 
of the local galaxy optical LF (eg iNorberg et al.ll2002h and 
the low mass e nd of the local mass f unction of supermassive 
black holes fee lShankar et al.ll2008bh . Therefore, if the radio 
LF extrapolated too far, then the calculated space density 



4 THE EVOLUTION OF THE MECHANICAL 
LF 

To go beyond the local Universe, we need to know the evolu- 
tion of the radio LF with redshift (assuming that there is no 
redshift depen dence in E qs. [2]and[ 4ll. We use the radio LFs 
determ ined bv lDunlop fc Pcacoclj (|l990l ) and lWillott et ahl 



Dunlop fc Peacockl f\199(i ) modelled the 2.7 GHz LF 



with the sum of three contributions: steep-spectrum radio 
sources in early-type galaxies, fiat-spectrum radio sources 
in early- type galaxies, and radio sources in late- type galax- 
ies. The latter are mainly powered by star formation and are 
thus irrelevant for our analysis. Flat-spectrum radio sources 
are beamed (the jets are aligned with the line of sight); 
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Figure 1. The thick solid lines show the bolometric LF of AGN inferred by I Hopkins et al from X-ray, o ptical and infrared data. 

The thin solid and dashed lines show the mechanical LFs inferred from the radio LFs of iDunlop fc Peacockl l ll99d) and IWillott et al] 
ll200lh . respectively. At ^ ~ 0, we have also shown the mechanical LFs inferred from lBest et al. I ll2005d) 's local radio LF (dotted lines). 
For each radio LF, we have determined two mechanical LFs, according to the two different fjct(^radio) conversions derived in Section|2l 
the red lines arc for Pjot oc L^^^^ (Eq. |4ll and the blue lines arc for Pjct oc L^^^^ (Eq.[2|. 



© 2005 RAS, MNRAS 000, [iHZ] 



On the jet contribution to the AGN cosmic energy budget 5 



1.00 



0.10 



E 
Q. 



0.01 



Synchro 


z=0 : 






-lO^Q ^10 ^10 '^lO ^10 ^10 "^10 

Lr.i„/L. 



10- 



Q. 



10' 



10^ 



E 



10- 



■ nOpKlnS eT ui. zuu/ 




■ 




7-'" 


Synchro 











Figure 2. The ratio between the mechanical power of radio 
sources (Eq. |6ll and the radiative power of luminous AGN (Eq. [7} 
in a representative volume of the local Universe. The solid line is 
computed using the mechanical power estimated from the size of 
the X-ray cavities (Eq. |4]|, in which case Pmoch depends on the 
minimum luminosity Lmin which one can extrapolate the LP in 
Eq. l(5j . The dashed line is computed using the mechanical power 
estimated from the minimum energy that is needed to explain the 
observed synchrotron emission (Eq. [2]l , in which case Pmcch docs 
not depend on this uncertainty. B oth lines are computed using 
the local LF of lBest et al. I l l2005al '). The vertical dotted line is a 
conservative lower limit for L„^]„. 



therefore, the equations in Section 2 overestimate their me- 
chanical luminosity. The addition of the LF of flat-spectrum 
radio sources to that of steep-spectrum radio sources makes 
httle difference to the latter, and so it is safe to deal with 
this complication by considering only steep-spectrum radio 
sources. The 2.7 GIfz LF of steep -spectrum radio sources in 
early-type galaxies computed by lDunlop fc Peacockl (|l990l ) 
is a double power-law function of the form of Eq. ([5]). In 
their pure luminosity evolution model (their other models 
give similar LFs out to 2: ~ 2) the dependence of redshift is 
entirely contained in L, = 10'^^-^^+'^-^'^''-° '^'^'^ W Uz'^^sr-^ 
We take this LF, convert it into a 1.4 GHz LF by assuming 
an Of = —0.8 spectral index, and correct for a cosmology with 

= 0.3, Ha = 0.7 and h — 0.7 so that we can compare 
our results with th e bolometric LF determin ed by Hopkins 
et al. (2007), since iDunlop fc Peacockl (|l990D had assumed 
Jljv/ — 1, Qa ~ and h — 0.5. The blue and red solid 
lines in Fig. 1 are obtained by convolving the 1.4 GHz LF 
determined in this manner with Eqs. [5] and |3J respectively. 
Both the blue and the red curves are significantly below the 
bolometric LF estimated by Hopkins et al. (2007; thick solid 
lines) at all redshifts. 

To check this result with an independe n t dete rmination 
of the radio LF, we consider Iwillott et all (|200lf l's best fit 
to the 151 MHz LF (their model C). We transform it into a 
1.4 GHz LF by assuming a = —0.8 and correct for the cos- 
mology (Willot et al. 2001 had assumed the same cosmology 
as Dunlop & Peacock 1990). The blue and the red dashed 
fines in Fig. 1 are obtained by convolving the 1.4 GHz LF 
determined in this manner with Eqs. [2] and |4j respectively. 
The dashed lines and the thin solid lines of the same colour 
run very close to each other. This demonstrates that the ra- 
dio LF is not a major source of uncertainty when it comes 



Figure 3. The mechanical (lines) and radiative (symbols) power 
of AGN per unit cosmic volume as a function of redshift. The 
mechanical power per unit volume p'^^^^^^° (dashed line) is com- 
pu ted by inserting Eg . ll2t i nto Eq. ^ and using the radio LF 
of iDunlop fc Peacockl lll99Clh . The solid line shows the redshift 
dependence of pj^gj.jj. The vertical normalisation of this line de- 
pends upon the lower luminosity limit of the radio LF adopted for 
the integration interval in Eq. but the redshift dependence 
does not. The plotted solid line uses a normalisation that cor- 
responds to Ljnin = 10^^ WHz~^ at z = 0. The radiative power 
per unit v olume Pi-^d is compute d from Eq. (0 using the bolomet- 
ric LF of iHopkins et al] l l2007t) . The dotted lines correspond to 

^'''''mcch"" ^'^'^ ^^PInach^ '^^'^V ^^'^ added to ease the comparison 
of the dashed line and the solid line with the symbols. 



to determining the mechanical LF of AGN, even at high 
redshifts. 

The mechanical power per unit volume obtain e d inte - 
grating Eq. ([2} over the LF of iDunlop fc Peacockl (|l990t ). 
Pmcch^°' shown as a function of redshift in Fig. 3 (dashed 
line) and compared to the radiative power per unit volume 
from Eq. (O (symbols). The dotted line in Fig. 3 is simply 
'^'^ Pmcch'° ■ shows that the cosmic evolution of the me- 
chanical power density traces that of the radiative power 
density, at least for Eq. ([2]). It also shows that throughout 
cosmic time jets contribute to a small fraction (~ 2%) of 
the AGN energy budget. This fraction goes up by a fac- 
tor of a few at low redshifts, if we believe that Eq. Q is a 
more accurate determination of Lmcch (solid line in Fig. 3). 
Even in that case, it is unlikely that the mechanical energy 
accounts for much more than ~ 10% of the AGN cosmic en- 
ergy budget locally. In addition, the cosmic evolution of the 
jet mechanical energy estimated by Eq. (4) is weaker than 
that of the radiative power, and so in this case the fraction 
of the AGN energy budget associated with jet mechanical 
power falls with increasing redshift, being well below 10% 
at 2 > 1. 



5 DISCUSSION AND CONCLUSION 

Fig. 1 suggests that the radio LF is not a major source of un- 
certainty when it comes to determining the mechanical LF 
of radio sources. It should be cautioned that constraints on 
the evolution of the faint end of the radio LF beyond 2 ~ 1 
remain quite poor. The uncertainties are almost certainly 
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larger than the variations between the dilTerent determina- 
tions of the radio LF evolution suggest. Nevertheless, it is 
also clear that the main source of uncertainty is the expo- 
nent of the Lmoch oc ifjjjjio relation. 

The argument based on X-ray cavities suggests that the 
importance of the mechanical-energy output decreases in lu- 
minous high-accretion-rate objects. On the contrary, the ar- 
gument based on the minimum energy to produce the ob- 
served radio emission suggests that the mechanical-energy 
output traces the luminous output. It is not s urprising that 
the la tter conclusion follows from Eq. ([2} since lWillott et ahl 
(|200ll '). from which Eqs. HJ and (0 are derived, find that 
the jet mechanical luminosity is proportional to the optical 
narrow-line luminosity. We speculate that the results ob- 
tained with the two methods may be different because the 
jets in powerful radio sources advance supersonically and 
often pierce through the ambient hot gas. In that case, the 
cavities are not inflated gently, therefore the work done to in- 
flate the cavities is much larger than pV. Important progress 
would come from observational studies comparing the me- 
chanical luminosities derived with the two methods for the 
same objects. 

Independently of whether one favours Eq. ([2]) or Eq. (|4]), 
Fig. 1 shows clearly that the bolometric luminosity func- 
tion is larger than the mechanical luminosity function by 
at least one order of magnitude at all luminosities. Eq. ((2} 
establishes a ~ 2% value at all redshifts for the mechani- 
cal contribution to the AGN cosmic energy budget (Fig. 3). 
Eq. Q gives a higher value, which depends on the minimum 
luminosity at which the radio LF levels off (Fig. 2). While 
this luminosity is uncertain, we can reasonably estimate that 
the contribution inferred from Eq. Q is larger than the con- 
tribution inferred from Eq. ([21 by a factor of ~ 3 — 10. It 
is thus unlikely that mechanical energy accounts for much 
more than 10% of the AGN cosmic energy budget with 20% 
as a firm upp er limit. This is in broad agreement with pre - 
vious studies (| Merloni fc Heinj|2008l : IShankar et aLll2008al ) , 
which is encouraging given the different methods that we 
have used in our analyses. 

This result implies that radiatively inefficient accretion 
is unlikely to contribute to much more than 10% of the BH 
mass in the Universe, unless the overall energy efficiency of 
accretion in this mode is substantially lower than the en- 
ergy efficiency in the radiatively efficient mode, ie unless a 
substantial fraction of the energy is advected onto the BH, 
rather than coming out as either photons or je ts. Based on an 
ADAF (advec tion dominated ac c retion flow; iNaravan fc Yil 
I 1994 ) model, iMerloni fc Heinj (|2008l ) suggested a kinetic 
efficiency for the production of jets of ek — 0.005 (where 
imech = tkM»(?) in low-accretion rate AGN, and thus con- 
cluded that ~ 18 — 2 7% of the BH growth oc curs in a radio- 
jet-producing mode. IShankar et al.l (|2008al 'l found a simi- 
lar value (~ 20 — 30%) by deriving a kinetic efficiency of 
Efe ^ 0.01 for the production of radio jets in radiatively ef- 
ficient radio-loud AGN, and by adopting this value for all 
radio sources. However, it is not clear whether such a low 
value for is appropriate for radiatively inefficient AGN, 
which may channel most of t he accretion power into jets 
iBlandford fc Begelmanll 19991 ). 

How do our results fit in the emergent scenario, in which 
jet heating plays a major role in the evolution of early-type 
galaxies and galaxy clusters? The energy released by the for- 



mation of a supermassive BHs is two orders of magnitude 
larger than the host galaxy's binding energy, so the issue is 
not the energy, but the efficiency with which it can absorbed 
by the ambient gas. Photoionisation of the inner orbitals of 
metals and Compton scattering are the main processes by 
which AGN radiation heats the surrounding gas. In nearby 
massive elliptical galaxies such as the systems studied by 
I Allen et al.l (|2006h , the gas on a galactic and group or cluster 
scale is hot and highly transparent, having uh ~ 10^^ cm~^. 
This column density is only a thousandth the column den- 
sity above which the gas becomes Thomson thick (ctt 
is Thomson cross section for electron scattering), meaning 
that only ~ 1 photon in 10^ is scattered by an electron be- 
fore leaving the galaxy. Even discounting any inefficiencies 
in the transfer of the photon energy to the gas in the scat- 
tering process (ie Compton scattering transfers ^'"^^i < 1% 
of the photon energy to a free electron per scattering event), 
the fraction of the luminous energy that would be absorbed 
by the gas if a quasar switched on in a nearby massive el- 
liptical galaxy would be ^ 0.1%. Therefore, based on this 
argument, mechanical heating could be > 10 times more im- 
portant than radiative heating even if the mechanical power 
were 100 times smaller than the radiative power. 

There is broad observational evidence that mechani- 
cal heating by jets plays an important role in solving the 
cooling-flow problem in galaxy clusters. The evidence that 
the same solution can be applied to individual galaxies is 
much weaker because jets may be coUimated on kilopar- 
sec scales and transport most of the energy to beyond the 
gaseous halo of the host galaxy. However, the bulk of the 
jet mechanical power is produced in low luminosity radio 
sources, most of which have small radio sizes. Even in larger 
sources, jet-interstellar medium interactions may also occur 
on sub-kiloparsec scales (the knots in the jets of M87 may 
be evidence for this). Therefore, in individual ellipticals the 
role of weak radio sou rces vs. episodic quasar heating (eg 
ICiotti fc Ostrikeill2007l 'l remains an open problem. 
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